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My early papers [7,8] on what is now called computed tomography were en- 
titled "Representation of a Function by its Line Integrals, with some Radio- 
logical Applications”, and in one of them I wrote "One would think that this 
problem would be a standard part of the nineteenth century mathematical reper- 
toire, but the author has found no reference to it in standard works". I have 
still found no nineteenth century reference to the problem, but in the years 
since 1971 I have come across a number of early discoveries of Radon's problem 
in the physical sciences, which I would like to share with you since they are 
not generally known. For more mathematical references, see John [21] and 
Helgason [19]. In addition, I will sketch my connection with the problem and 
give some results which I have recently obtained. 

The first person I know of who tackled Radon's problem was the great 
Dutch physicist H.A. Lorentz. He found a solution to the three dimensional 
problem where a function is to be recovered from its integrals over planes. 

If £(p, 1) is the integral of f over a plane perpendicular to the vector pn from 
the origin, and a distance p from the origin, then f at the origin is given by 
g > 

£(0) = - oe au 
where de, is an element of solid angle in the direction of ne Since the origin 
may be arbitrarily chosen, the result holds for any point. We have no idea : 
why Lorentz thought of the problem, or what his method of proof was, and the 
only reason we know of his work is that the above result was attributed to him 
by Bockwinkel [4] who used the result in a long paper on the propagation of 
light in biaxial crystals. 

Following Lorentz came Radon's paper, but since all of you have copies of 


it, I shall not comment on it. 
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The problem appears again in Holland in 1925 and again involves the 
physicists Paul Ehrenfest and George Uhlenbeck (who is still alive in this 
country). In Uhlenbeck's paper he says that Ehrenfest drew his attention to 
the result of Lorentz and suggested that he generalize it to n-dimensions, 
which he proceeded to do using Fourier techniques [33]. No reason for solving 
the problem was given, and my guess is that it was just a cute problem crying 
out for a solution. I am indebted to Alberto Grunbaum for drawing my attention 
to the Bockwinkel and Uhlenback papers. 

1936 was a vintage year for it contained two independent discoveries of 
Radon's problem. One was in Stockholm where Cramér and Wold [15] asked the 
following question: given all marginal distributions of a probability distri- 
bution, can one infer the distribution itself? Marginal distributions are what 
people in computed tomography call projections or views, so the question is 
just Radon's problem, and in solving it, Cramér and Wold again used the Fourier 
integral approach. The second was in Leningrad where the well~known Armenian 
astronomer, V.A. Ambartsumian, posed the following question which he attributes 
to Eddington [3]. If one looks in a particular direction in space, one sees 
many stars and these are moving about relative to one another and to the sun, 
Astronomers would like to know their distribution of velocities but, except for 
a few close stars, one can only measure their radial velocities, which are de- 
duced from the Doppler shifts cf their spectra. The problem then is to deduce 
the actual distribution of velocities in three dimensions in space from the 
distributions of radial velocities in various directions. (It is, of course, 
necessary to assume the velocity distribution to be invariant under spatial 
translations.) This is just Radon's problem in a three dimensional velocity 
space rather than ordinary space, and Ambartsumian gave the solution in two and 
three dimensions [1] in the same form as Radon. Furthermore, he took groups of 
stars of three different spectral types, with four to five hundred stars in 
each group, and he used his theoretical results to deduce their actual velocity 
distributions from the distributions of their radial velocities [2]. One of 
his results is shown in Fig. 1 and it is in fact the velocity distribution pro- 
jected onto the galactic plane. This is the first numerical inversion of the 
Radon transform and it gives the lie to the often made statement that computed 
tomography would be impossible without computers. Details for the calculation 
are given in Ambartsumian's paper, and they suggest that even in 1936 computed 
tomography might have been able to make significant contributions to, say, the 
diagnosis of tumors in the head. 

The problems of finding abnormalities in the soft tissues from X-rays, 
even with the aid of computed tomography, are far greater than most people 


realize so it seems to me quite possible that Ambartsumian's numerical methods 
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Velocity distribution of stars projected onto the galactic plane. 


might have made significant contributions to that part of medicine had they 
been applied in 1936. I am indebted to Professor Ambartsumian for drawing my 
attention to his work and for telling me that he was informed of Radon's 
results two years after he published his work. 

In 1947 there is a fascinating abstract by Szarski and Wazewski [32] 
which begins: "En disposant des fotographies roentgeniques de la tête humaine 
de tous les côtés, peut-on déterminer la structure intérieure du crane, voici 
un problème proposé par M. Stanislaw Majerek, médicin de Cracovie qui, sous 
un forme élémentaire, a indiqué un méthode de solution qui nécessite une justi- 
fication mathématique". They then go on to describe Radon's problem and formu- 
late it in terms of a set of “fonctions cylindrique" and state that the problem 
consists of finding whether this set of functions tends to a solution of Radon's 
problem. I have been unable to find any further references to this work. 

If one looks at the sun, as Bracewell was doing around 1956, with a long 
radio antenna with a parabolic section perpendicular to its length the signal 
one gets is an average over a thin strip of the sun of its radio emission at 


some wavelength. On moving the antenna one gets averages over different strips, 
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and once again we have Radon's problem, which Bracewell proceeded to solve 
using Fourier and other methods [5]. 

As I have recounted elsewhere [10] in 1956 I came across the problem in 
radiology where one has to determine the X-ray absorption coefficient from 
point to point in a disc from attentuation measurements on X-ray beams travers- 
ing the disc; Radon's problem again! Since the treatment I gave differs from 
previous treatments, and since it shows some interesting features which appear 
later, I will outline it. 

First consider a circularly symmetrical function f(r). Then if £(p) is 
the line integral along a line which is a perpendicular distance p from the 
origin, it is easy to show that 


a f(r)rdr 
f(p) = 2 | ===. (1) 
prep 
This is just Abel's equation with the well-known solution 
1 £' (p)dp 
f(r) = - = | — (2) 
T la PA 


which is quite amenable to numerical analysis and which I used in an experiment 
[7]. The fact that the integral in Eq. (2) extends from r to © is an explicit 
expression of the hole theorem: to find f(r) one needs only Ep} for p > r. 
This remains true for non-circularly symmetrical functions as can be seen from 
the following. Let the line L along which f is calculated be specified by 
(p,ġ) as in Fig. 2. 


Reconstruction of Densities from their Projections 


Fig. 2. Coordinate system defining the line L in u slice of the object. The line L is the 
limit of the cylinder C of fig. 1. 


Figure 2. 


Coordinates for the Radon transform in R° 


Then if one writes 


£(r,0) = E f(e, 80,9) = 1 £,(pde 
2 2 


ikp 
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one can show that 


m © £,(r)T, (p/r) rdr 
f) = f r , (3) 
P 
which has the solution A 
f; (p)T, (p/r)dp 
goat 2s 
f(r) =- 5 [ ner, , (4) 


where Ty (x) is the Tschebicheff polynomial of the first kind which is orthogonal 
in (-1,1) [8]. In Eq. (4), T, (x) appears for x > 1, and its leading term is 
const X e hence if fy (P) is obtained from noisy data, the noise is propagated 
badly from the exterior in as one goes to smaller values of r. In a recent 
application of this method in fusion physics [30] it is shown that Eq. (4) 

gives bad results for & R, 10. Hence Eq. (4) is not good for numerical solu- 
tions in medicine where one would require 2 to be of the order of hundreds. 

I therefore resorted to radial expansions of fo and A which take the following 
form. If 


E (E) = aa) o (A, t, r°), (5a) 
then 
fo (p) = const * (1 29112 Ceo im. (5b) 


where the Gk are Jacobi polynomials in (0,1) and the ch are Gegenbauer poly- 
nomials. This expansion, for general A, was recently given by Louis [22] and 
I used only the case A = 1 which has been discussed by Marr [24], Ein-Gal [17] 
and others. Application of these formulae in an experiment to find f from f 
showed that computed tomography would work well for real data, but today the 
most commonly used numerical methods are based on the convolution algorithm of 
Ramachandran and Laksminarayanan [29]. 

When Radon's problem is extended to R° Ludwig [23] and Deans [16] 
obtained formulae which contain Eqs. (3) and (4) as special cases. For more 
general spaces, a more general form of the hole theorem is immediately obtain- 
able from a theorem (2.1) of Helgason [18]. 

Since the middle sixties there has been an explosion in applications of 
Radon's problem to many problems in science and engineering some of which are 
given in Marr [25] and more in a hard-to-obtain bibliography by Richard Gordon. 
Since the advent of commercial computed tomography scanners [20] there has been 
a greater explosion in the medical field, and reviews of Brooks and DiChiro 
[6] and Shepp and Kruskal [31] are useful, though both are becoming dated. 

Radon's problem in R? gives rise to a nice tiling problem which has not 
been used and which may not even be useful. Consider two pairs of points, 
S]»5, and s5 which define a beam of X-rays traversing a sample in the follow- 


ing way: only lines which intersect S153 and S153 contributed to a measurement 
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of some average value of f with this beam. The lines defined by S,S, and ss 


are represented by a region of the Radon transform bounded by four círcles, 


2 


namely, those with radii 05,» 05,» Osi, 05, where O is the origin of the 
transform. The question then is how the Radon transform of a finite object can 
be covered exactly by a set of such beams which intersect it. The solution is 
easy for a model of positron-annihilation (PET) scanning when the detectors can 
be considered to be arcs of the circle defining the object [11]. For exactly 
parallel beams (i.e. S,S, and S¡S, forming the corners of a fixed rectangle) 
there is no exact solution, but a mix of parallel and almost-parallel beams 
does provide a solutíon [9]. One wonders whether there are other solutions, 
and whether these results can be extended to R”. 

When I had obtained a solution of Radon's problem in the form given in 
Eq. (4), I asked myself whether there were other curves specified by parameters 
(p,¢) for which Radon's problem could be solved. I immediately tried circles 
through the origin with diameter p and orientation of the diameter relative to 


the x-axis given by $, and obtained results remarkably similar to Equations (3) 


and (4), namely 
p £,(r)T, (r/p)dr 


fo (P) = 2 A (6) 
which has the solution 


r f,(p)T,(r/p)dp 
Pp 2 . m 


E -5 = r°-p 
Quinto and I recently generalized these to R" and applied the results to solu- 
tions of Darboux' equation (14]. The results will not be given here, but the 
kernels are more complicated than those in Eqs. (6) and (7), involving Cegen- 
bauer polynomials which appear here as in [23] and [16] as a result of using 
the Funk-Hecke theorem for spherical harmonic expansions. 

As an older physicist, I was aware that in electrostatics straight lines 
and planes are related to circles and spheres by inversion ((r,0) > (1/r,0)) 
and this connection intrigued me for a long time. Recently I came across two 
families of curves which are related by inversion, and for which Radon's problem 
can be solved, which include the straight lines and circles. These are, in 


the previous notation, 


r* cos{a(6-¢)} =p", a> 0 
(8a) 


A, B > 0, 


e cosí8(8-4)) 


For a = 1/2, 1, 2 we have parabolae, straight lines and hyperbolae respectively, 
and for ß = 1/2, 1, 2 we have cardioids, circles through the origin and one 
branched lemniscates of Bernoulli, respectively. Many of these curves are the 


orbits of particles in central fields or of charged particles in certain 
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magnetic fields. The details can be found in [12], but in general the results 
are very like those for the ordinary Radon transform: the inversion formulae 
are similar and the appropriate hole theorems make their appearance. More 
recently [13] certain other properties have been obtained. For general values 
of a and B certain inverse and direct powers of r (r # 0) yield null Radon 
transforms, as pointed out by Newman [26] and Louis [22] for the ordinary 
Radon transform. For a = 1, 8 = 1, Quinto [28] has shown that these constitute 
the entire null space. For a or 8 = 1/m, m= 1, 2, 3... a number of orthogonal 
expansions of fo and E, are given. For a-curves in (0,1) the expansions are an 
extension of the results of Louis (Eqs. 5a,5b); in (1,%) they are an extension 
of Perry's results [27]. In (0,%) they follow the pattern given in [8]. For 
B-curves, the inversion property allows one to write down immediately the 
appropriate formulae for (0,1) and (1,®) from the results for a-curves for 
(1,©) and (0,1) respectively. 

Working on Radon's problem has given me a great deal of pleasure and, as 
usual, each question answered has raised new questions. There are still plenty 


of these to be answered. 
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